Thermopower of correlated semiconductors : application to FeAs2 and FeSb2 
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We investigate the effect of electronic correlations onto the thermoelectricity of semi-conductors 
and insulators. Appealing to model considerations, we study various many-body renormalizations 
that enter the thermoelectric response. We find that, contrary to the case of correlated metals, 
correlation effects do not per se enhance the Seebeck coefficient or the figure of merit, for the former 
of which we give an upper bound in the limit of vanishing vertex corrections. For two materials of 
current interest, FeAs2 and FeSb2, we compute the electronic structure and thermopower. We find 
FeAs2 to be well described within density functional theory, and the therefrom deduced Seebeck 
coefficient to be in quantitative agreement with experiment. The capturing of the insulating ground 
state of FeSb2, however, requires the inclusion of many-body effects, in which we succeed by applying 
the GW approximation. Yet, while we get qualitative agreement for the thermopower of FeSb2 at 
intermediate temperatures, the tremendously large Seebeck coefficient at low temperatures is found 
to violate our upper bound, suggesting the presence of decisive (e.g. phonon mediated) vertex 
corrections. 
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I. INTRODUCTION 

While recent efforts to design materials with enhanced 
thermoelectric properties were mainly focused on reduc- 
ing the lattice contributions to the thermal conductivity 
by super- or nano structures,— interest in the potential 
merits of electronic correlation effects was revived by the 
discovery of large Seebeck coefficients in transition metal 
compounds, such as FeSi;^ Na:rCo02(^ and FeSb2^, with 
the latter displaying an astonishing response of up to 
S = -45 mV/K at 12Ki. Indeed, on a model level^"^ 
as well as for realistic compoundS)2ii£ correlation effects 
were shown to enhance the Seebeck coefficient in metals 
and transition metal oxide o^^i*^ . 

In this work we address the thermoelectric response of 
correlated semiconductors and insulators. Our aim is to 
set up the general formalism for discussing the thermo- 
electric response, with the ultimate goal to understand 
the origin of the large thermoelectricity observed in corre- 
lated semiconductors and to search for high performance 
thermoelectrics in this class of materials. This work par- 
allels the analysis done for correlated metals in Refsi^ii^. 

In particular, we investigate whether the electronic 
structure and correlation effects alone can account for 
the very different magnitudes in the Seebeck coefficient 
of the two iso-structural and iso-electronic compounds 
FeSb2 and FeAs2. In the case of FeSb2, it has indeed 
been conjectured that electronic correlations are at the 
origin of the huge thermoelectric response^ii^"— . 

The setup of the paper is the following. First we de- 
scribe the FeX2 materials. Then we shall first extent gen- 
eral text-book considerations for the Seebeck coefficient 
to include the important aspect of carrier selective elec- 
tronic renormalizations, as well as properties beyond the 
picture of coherent band-structures. We use this frame- 



work to make general arguments on how to obtain high 
values for the figure of merit ZT. Given the sizes of 
the charge gaps, effective masses, and other parameters, 
these considerations allow us to put constraints on the 
possible regimes of materials of interest. In Section lTVl we 
apply realistic electronic structure tools, and show that 
for iron arsenide, FeAs2, both the electronic structure, 
as well as the thermoelectric response can be understood 
as that of a conventional semiconductor and described 
quantitatively by ah initio band-structure methods. For 
the iron antimonide on the other hand, standard density 
functional theory (DFT)i^ based methods are known to 
be insufficient to account for the electronic structur e* ''i^^ . 
To show that correlation effects play an important role, 
we employ a hybrid functional approach^^, as well as 
Hedin's GW approximation^^, with the latter yielding 
results in good agreement with the experimental charge 
gap. As to the thermoelectric response, however, we find 
qualitative agreement with experiment only at intermedi- 
ate temperatures (35-70K). Our analysis shows that the 
low temperature Seebeck coefficient of FeSb2 is incom- 
patible with a local electronic picture, suggesting the im- 
portance of vertex corrections and non-local self-energy 
effects (that we neglect), or the presence of a substan- 
tial phonon drag effect^ as is e.g. found in the classical 
example of p-type Germanium^^. 



II. THE MATERIALS 

Despite the structural similarity of FeSb2, and FeAs2, 
experimental findings point to markedly different proper- 
ties, heralding a varying importance of correlation, and, 
potentially, electron-phonon effects. 

FeAs2 is an insulator with a gap of 0.2-0.22 eV— i^, 
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as obtained from the activation behavior in the resistiv- 
ity at temperatures of 200K and higher. Below 200K, 
the influence of impurities is pivotaU^^ : the resistiv- 
ity has a metallic slope before resuming, below 30K and 
down to lOK, an activation law with an energy of O.OleV. 
Further, below lOK, the resistivity exhibits activation 
with 6K (O.SmeV).— The Hall coefficient is negative for 
all temperatures^^. Congrucntly, the Secbcck coefficient 
is negative as well : From its room temperature value 
— 200/xV/K^. it grows in magnitude upon cooling, to 
reach — 7mV/K at 12K,— before it vanishes towards zero 
temperatureF^^ 

In the case of FcSb2 , optical spectroscopy finds a small 
gap of 432K (37meV) at low temperatures,— but wit- 
nesses the development of a Drude-like peak at 70K and 
above. The concomitant transfer of spectral weight is 
found to extend over an energy range as high as leV— , 
i.e. a scale that is much larger than the initial gap, a 
common harbinger of correlation effecta^^. The resis- 
tivity of FeSb2, on the other hand, has three distinct 
temperature regimes that exhibit activated behavior : 
In the range of 50-lOOK the activation energy corre- 
sponds to a gap of BOOK [2&uieY)^. From 20K down 
to lOK, the resistivity shows a shoulder-like behavior 
with an activation energy of A/2 = SmeV, while be- 
low 5K, extrinsic impurities are believed to be at the 
origin of a weakly temperature dependent resistivity fol- 
lowing an activation behavior with 0.04-0.09meV4. The 
resistivity is anisotropic4^i^, and some experiments find 
metallic transport behavior [dp/dT > 0) for selected 
directions above AOK^^^ This anisotropy is also seen 
in Hall measurements As a matter of fact, the Hall 
coefficient even changes sign for some polarizations (at 
lOOK,— or 40K— ), with predominant electron character 
(n-type) below these temperatures. The Seebeck coef- 
ficient at 300K is found to be 15 - 40/zV/K, with the 
sign depending on the polarization.i^^ Upon on low- 
ering the temperature, the Seebeck coefficient passes a 
local maximum {d^S/dT^ > 0) at around 40K, before 
turning towards very large negative values, reaching, de- 
pending on the polarization and the sample, a global ex- 
trcnia of up to — 45mV/K at slightly above lOKi^ Below 
this temperature, the coefficient drops sharply in mag- 
nitude and practically vanishes at 5K and below. In- 
terestingly, the largest thermopower is thus found in the 
temperature range where the resistivity has the shoulder- 
like behavior. Noteworthy, this regime is concomitant 
with the appearance of a prominent feature in the "elec- 
tronic" specific healji^. No such feature is found for the 
arsenide^. However, the unlocking of spins in FeSb2 
becomes appreciable only beyond this regime at around 
150K, where the entropy reaches i?log2, owing to a sec- 
ond and larger hump in the specific heat, and in con- 
gruity with the susceptibilit y^^ i . Indeed FeSb2 becomes 
paramagnetic above lOOK— , and a Curie-like downturn 
appears at temperatures above 350K— , whereas the sus- 
ceptibility of FeAs2 is flat up to 350Kii. That the low 
temperature feature in the specific heat of FeSb2 has no 



spin signature might indicate that its contribution to the 
entropy is associated with either the charge degrees of 
freedom or an clcctron-phonon effect. The importance of 
electron-electron effects in FeSb2 is further highlighted 
by the fact that various properties are very sensitive 
with respect to changes in the carrier density. Doping 
the system with electrons, e.g. FeSb2-i:Snaf^, or holes, 
e.g. FeSb2-a:Te2^i^, instantly metalizes the compound, 
sometimes generates a Curie law at low temperatures^ 
and reduces the Seebeck coefficien t^^i'^^ . Thermoelectric 
properties of FeSb2_a;Tea; are indeed that of a correlated 
metal, i.e. the low temperature Seebeck coefficient is lin- 
ear in T, with an enhancement factor of 15 via the effec- 
tive massif. 

Despite these indications for correlation effects, some 
experimental findings for the antimonide are quantita- 
tively reproducible by conventional band-structure meth- 
ods : Volume and bulk-modulus^ are very well captured 
within the generalized gradient approximation (GGA) of 
DFT— . Also, the finding of small electron and hole pock- 
ets in FeSb a^^i^^ (see also below) within band-structure 
methods could simply be attributed to the well doc- 
umented underestimation of charge gaps within DFT. 
Moreover, a calculation within the local density approx- 
imation (LDA) with a Hartree-like Coulomb interac- 
tion (LDA+U) suggested that FeSb2, while being para- 
magnetic, could be close to a ferromagnetic instabil- 
itjsii. Weak ferromagnetism was then indeed found in 
Fci_,Co,Sb2^. 



III. TRANSPORT FORMLUAE, GENERAL & 
MODEL CONSIDERATIONS 

Within the Kubo formalism the Seebeck coefficient - 
that relates the gradients of te mper ature and electrical 
field - is given by (see e.g. Ref. lolllOl) 

|e| Ao 

where the current-current (current-heat current) correla- 
tion function Aq (Ai) is given by 

An = J dwiriio-^^r (^-^^Eii,) (2) 

Here, is the Fermi function, n is the Fermi level, and 
S is the transport kernel. If vertex corrections are ne- 
glected, the transport kernel can be expressed (in matrix 
notation) as 

Eico) = J2Tr[vik)A(k,Lo)vik)A{k,co)] (3) 
k 

with the Fermi velocity Wij(k) = (^kJ Vji/'jj.p, and 
the spectral function j4y (k,a;), where ■0]^^ is a complete 
set of one electron basis functions, such as Kohn-Sham 
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orbitals. Using transport coefficients Eq. ([2]), we can fur- 
ther express (see e.g. Ref. 0) the dc conductivity, the 
thermal conductivity, and the figure of merit as 



hV 



2TTk% f , A? 



ZT = 



(4) 
(5) 
(6) 



where V is the unit-cell volume, and the thermal lat- 
tice conductivity. 

The chemical potential, ^, is obtained by the require- 
ment of the charge neutrality, 



n — p — ~ 



(7) 



where 



E 



is the number of electrons (holes). Here, A'^''"{k,uj) is 
the valence/conduction spectral function. We also al- 
lowed for the presence of ionized donor impurities njj+ = 
no [l + 2e~'^''^"~''^] of concentration no at an energy 
Ed- 



A. model considerations 

Here we find it instructive to extend on the usual text- 
book considerations (see e.g. Ref. [36l ) and generalize to 
include carrier dependent masses, renormalizations, as 
well as finite (yet energy-independent) scattering ampli- 
tudes. This makes it possible to investigate the impor- 
tant effects of particlc-holc asymmetry, carrier coherence, 
and allows us later to discuss the consistency of a purely 
diffusive thermopower for a given material. 

Assuming a Lorentzian line shape of the conduction 
(c) and valence (v) spectral functions 



(..-CT + (rr)' 



(9) 



i.e. we limit the discussion to the coherent part, Acoh, (of 
weight Zfc) of the full spectrum, A = Acoh + Aincohi with 
a quasi-particle dispersion ^fc, and an elastic scattering 
of amplitude Ffe. Within this approximation, one finds 
for the number of electrons^ 



^ + |^(r. + z(a-A^) 



(10) 



where tjj^z) is the digamma function, and /3 = l/{kBT) 
the inverse temperature. The above expression reduces 
to the usual J^k /(a* ~ Cfc) in the coherent limit (F = 0). 



The response functions, Eq. ([2]), can be also be expressed 
analytically (here, we restrict ourselves to n = 0, 1) as 

Y.^^W-i'TR^i^'k^) (11) 

i. ^ V 



An — 
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-|^((6 - M)"rfci?eV^"(x) nVllmr{x, 

where the arguments of the derivatives of the digamma 
function are 

i + ^(r.-.(a.-.) 

The first contribution in the curly brackets of Eq. \W\ 
is the leading term in the coherent limit (F — > 0). Indeed 
with /3/(27r2)i?eV''(l/2-|-i/3a;/(27r) = -f{x) one recovers 
the Boltzmann expressions for transport coefficients)^^ 
in which appears an ad hoc lifetime r = 1/(2F). In the 
above expression, however, the influence of finite scat- 
terings is not restricted to the pre-factor, but since the 
spectrum broadens, a wider energy range becomes acti- 
vated for supplying charge carriers, and, as a result, the 
Fermi statistics assumes the digamma form, and, also, 
higher order terms appear. Therewith, contrary to the 
Boltzmann description at small F, coherence effects do 
not in general cancel in the ratio Ai/Aq of the Seebeck 
coefficient. 

For illustrative purposes, we now analyze Eq. (jlip in 
terms of special cases for simple quadratic dispersions : 
We consider bare bands e^ .„ = , and an interacting 
dispersion 



= ±A/2± 



21,2 



2m* 



(12) 



where A is the charge gap, m* the effective mass of the 
carriers, and the origin of the chemical potential is chosen 
at the mid-gap point. Further, we assume the Fermi 
velocities to be given by the group velocity^ 

Vk = \dkek (13) 
n 

and weights, Z , and scattering rates, F, that are inde- 
pendent of momentum. 



1. the large gap cotierent semi-conductor 

In the limit of a coherent system (F <C 1) with a large 
gap (/3 [A/2 - > 1)7 Eq- CB can be simplified to 



At' = 

(±l)"^^e- 



(14) 



'^^/^|[/3(±/.-A/2)r-^n| 



where all carrier specific parameters have been gathered 
in 
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(15) 
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(a) coefficient of the 1/T behavior in the gap. (b) maximum value of the 1/T coefficient. 



FIG. 1: The symmetric, coherent, large gap semiconductor, (a) coefflcient of the 1/T low temperature behavior of the Seebeck 
coefficient as a function of the chemical potential, and for different temperatures, (b) temperature evolution of the extremal 
value of the 1/T coefficient of the thermopower. 



Therewith the Seebeck coefficient becomes 
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(16) 
(17) 



where the asymmetry parameter SX (that depends on fi 
and T) is given by 



6X 



- A" 
X^TX" 



(18) 



Hence, a large Seebeck coefiicient can be achieved by an 
interplay of the gap. A, the anisotropy or asymmetry 
6X in the transport function - stemming from either the 
densities of states (m^„), different bandwidth narrow- 
ings (to*^), scattering amplitudes (r^^) or quasiparticle 
weights {Zc,v)- 



2. Upper limit for thermopower in a semiconductor 

The position of the chemical potential plays an impor- 
tant role in maximizing the thermopower of a semicon- 
ductor. Let us take the example of two equivalent bands. 
In this case, the term in brackets in Eq. p7)) becomes 



^ - A/2tanh(^^) 



(19) 



This coefficient of the 1/T behavior is displayed in 
Fig.[TJa) as a function of the chemical potential for a few 
different temperatures. At high temperatures the opti- 
mal chemical potential, which maximizes thermopower, 
is near the gap edges. At low temperature ksT < A/ 2, 
the 1/T coefficient shows a local extrema. The Seebeck 
coefficient vanishes at the point of particle-hole symme- 
try (/^ = in the symmetric case considered here — ) , 
and the optimal location of the chemical potential at low 



temperature is in the direct vicinity of that point, hence 
very close to the center of the gap. 

Fig.[ljb) displays the value of the 1 /T coefficient in this 
extremum and the value at the gap edge as a function of 
temperature. For fc^T ^ 0.3A the maximum value of the 
thermopower is achieved when the chemical potential is 
at the gap edge, and for lower temperature, it is achieved 
close to the middle of the gap, where the thermopower 
can reach the maximum value of 5* = A/(2eT). 

In an asymmetric case, the thermopower can be larger 
than this maximum value, however, for a given charge 
gap A, there is always an upper bound for the Seebeck 
coefficient, namely 



\S{T)e\ < A/T + 5/2kB. 



(20) 



This is because the asymmetry is bounded to an absolute 
value of one \6X\ < 1. This extremal value corresponds 
to the fictitious system in which only one type of the two 
carriers contributes to the thermoelectricity, e.g. the con- 
duction electrons, and with the chemical potential being, 
in that case, at the edge of the valence band. 

Thus the correlation effects, such as small Z in the 
conduction band and large Z in the valence band, can 
enhance the thermopower of a semiconductor. However, 
this effect is limited by the form of Eq. P7)) allowing 
maximum S bounded by Eq. (j20p . The possible mer- 
its of electron-hole-asymmetry for the case of metals is 
discussed in Ref. llOl. 



3. Model semi-conductor in the presence of donor 
impurities 

With the goal of understanding the thermopower and 
the figure of merit in a renormalized semiconductor in 
a very general setup, including the presence of impu- 
rities, we now numerically study the model based on 
the response functions Eq. (HJ). As before, we assume 
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parabolic dispersions Eq. p2|) . with the band structure 
depicted in the inset of Fig.[5Jb) : excitations of different 
effective masses are separated in energy by a gap A, and 
we allow for the presence of donor impurities, situated at 
an energy Ed, as measured from the middle of the gap. 
We again assume transition matrix elements to be given 
by the group velocity, Eq. ([T5|) . 

We choose the parameters compatible with the band 
structure of FeAs2: we consider a gap A = 0.2eV, 
and, unless stated otherwise, an impurity level at Ejj = 
95meV, as inferred from the low temperature activation 
behavior of the resistivityi^. In our current treatment, we 
assume that the impurity carriers have vanishing Fermi 
velocities, and thus their only effect is to shift the chem- 
ical potential. 

To fix the particle-hole asymmetry, we note that at 
high temperatures, the number of ionized impurities are 
irrelevant with respect to the number of conduction and 
valence carriers, and the chemical potential follows the 
intrinsic behavior. In the coherent limit of the large gap 
semiconductor (see above) one finds that 

/i = 3A:BT/41n(?/„/77c), with t] ^ m* /mo- (21) 

In this regime the resistivity shows an activation law with 
the activation gap A/2. Of course, the situation in a real 
material can be much more complicated (several types 
of impurities, temperature dependence of the gap, etc.). 
Using our ab initio data for FcAs2 (that is presented be- 
low in Section ITVl and Section FV Ap . and assuming the 
approximate validity of Eq. (|2ip for non-parabolic disper- 
sions, one finds the ratio of the valence and conduction 
effective mass rju/ijc = 2.5 for FeAs2, which we will use 
for all the following model calculations. 

Further, we note that a uniform weight-factor Z can- 
cels in the Seebeck coefficient, whereas in the figure of 
merit it can be seen as a scaling factor of the ther- 
mal lattice conductivity, for which we assume kl jZ"^ = 
250W(Km)~^j^ Moreover, we use a unit-cell volume of 
80A3. 

Having thus fixed the size of the gap and the asymme- 
try, the principle parameters to vary in this setup are the 
concentration of impurities nu and the scattering rate 
r (that we assume to be orbital independent). We will 
also study the dependence on the position of the impurity 
band Ejj away from the value motivated by the experi- 
mental resistivity. 

a. impurity concentration. With these parameters, 
we display in Fig. [DJa) the temperature dependence of 
the chemical potential, the Seebeck coefficient, and the 
figure of merit for various impurity concentrations no, 
and for a constant scattering rate F = 5/ieV. 

In the intrinsic case, uq = 0, the chemical potential is 
indeed linear above a certain temperature that is related 
to the scattering rate F. Below this regime, the chemi- 
cal potential is almost temperature independent. Since 
r is small, the point of charge neutrality at zero temper- 
ature is very close to the midgap point (where it is in the 
coherent case). The point of particle- hole compensation 



for the thermopower, however, is higher in energy than 
the midgap point, as rationalized, above, for the sym- 
metric semiconductor. As a consequence, the chemical 
potential is actually moving through the symmetry point 
of the Seebeck coefficient, and we expect it to change 
sign as a function of temperature. This is clearly seen in 
Fig.[2Ia). 

At finite dopings, and at low temperatures, the donor 
level plays the role of the valence band, and the chem- 
ical potential (top panel) is between Eo and the con- 
duction band. In a coherent semiconductor the chemi- 
cal potential is /i « (A/2 — S/2) — kBT\n{n/n'jj), where 
6 = A/2 — Eo is the impurity activation energy, seen 
in the resistivity of this regime (p ~ e~^^^'^). At very 
low temperatures, /x is pinned to Eo- At intermediate 
temperatures the chemical potential smoothly connects 
with the intrinsic high temperature slope, as shown in 
Fig. m^a) . For large enough concentration of impurities 
(no > 10^^/cTO'^), the chemical potential can even go in- 
side the conduction band at some intermediate temper- 
ature, which can result in a shoulder, or even a metallic 
slope in the resistivity (cf. the mentioned transport mea- 
surements on FeSb2, and FeAsa^^). 

For finite impurity concentration {no > 0), the See- 
beck coefficient displayed in Fig. [U^a) may or may not be 
enhanced at a given temperature, depending on whether 
or not the additional carriers bring the chemical potential 
closer to its optimal value. As explained above for the 
case of two particle hole symmetric bands, at fixed tem- 
perature there exists a value of the chemical potential, 
which maximizes the Seebeck coefficient. In our asym- 
metric setup with rj^/rjc > 1, the optimum chemical po- 
tential is located above the mid gap point. In the limit 
of vanishing impurity concentration, the midgap remains 
the point of charge neutrality at zero temperature, while 
the optimal chemical potential is very near that point, 
hence the Seebeck coefficient is a very strong function of 
temperature in this limit, and can even change sign, as 
seen in the middle panel of Fig. [H^a). We note however, 
that the fundamental extremum, established by the size 
of the gap, is always respected as is evident in Fig.^a). 

The efficiency of the thermoelectrical material is de- 
termined by its figure of merit ZT, which we also plot 
in the lower panel of Fig. HJa). ZT can be greatly en- 
hanced by the presence of impurities, and its maximum 
is not necessarily in close vicinity of the thermopower 
maximum. Indeed the largest ZT for the current param- 
eters is achieved at about T = 220K for a concentration 
no = 2 • 10^^/cm^. The Seebeck coefficient at the this 
point is actually smaller than in the intrinsic limit. The 
position of the impurity with respect to the gap edge has 
a large effect on the optimal impurity density. While for 
our specific choice of parameters, A/2 — Eo ~ 5meV, 
the optimal density is no = 2 • 10^^/cm'^, we notice that 
for larger gaps and/or larger separation of the impurity 
level from the gap edge, the optimal impurity density can 
reach a value as large as ~' 10^°/cm'^i^ To elucidate the 
origin of an optimal density further, we show in Fig.[2Ib) 
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(a) effect of the impurity concentration n^i on fi, S, and ZT. 
(A=0.2eV, £;D=95meV, r=5/ieV) 
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(c) effect of the scattering amplitude F (top) and the position 
Ed of the impurity (bottom). (A=0.2eV, no=2 ■ lO^^/cm-^*). 
Inset: S{T) for various scattering rates. 
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(b) carrier density dependence of 5, a, and k. (T=220K, 
A=0.2eV, r=5/xeV). Inset: Sketch of the quasiparticle 
dispersion. 
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(d) ZT with and without lattice contributions (top), carrier 
concentration that maximizes S or ZT (bottom left), ZT as a 
function of (bottom right). (A=0.2eV, r=5^eV) 



FIG. 2: Thermoelectric properties of the model semi-conductor with parameters inspired from FeAs2. (a) from top to bottom : 
chemical potential, Seebeck coefficient, and figure of merit ZT as a function of temperature for different impurity concentrations, 
(b) ingredients to the figure of merit ZT = S'^aT/n at T = 2207^ as a function of carrier density. The inset shows the general 
setup of asymmetric valence and conduction dispersions, with an impurity level Ed in the gap A. (c) effect of the scattering 
amplitude V (top) and the impurity level position Ed (bottom) onto the figure of merit ZT, and the thermopower (inset), 
(d) top ; comparison of the figure of merit from (a) with the purely electronic figure of merit [kl = 0), bottom left : carrier 
density that maximizes \S\ and ZT (for the chemical potential in the upper half of the gap), bottom right : ZT as a function 
of the donor concentration for several temperatures, and a fixed impurity level Ed = 95meV. We assume a unit-cell volume 
V = 80A^, and a thermal lattice conductivity (scaled with the quasi-particle weight Z, see text.) of kl/Z^ = 250W/(Km). 



the dependence of the quantities entering the expression of ZT as a function of the particle density (n = p+n^) at 



7 



fixed T = 220K, the temperature which maximizes ZT in 
Fig.[2{a). We use here the total particle density, because 
in this case the description becomes independent of the 
bare impurity concentration and the level position, and 
in particular the additional carriers can find their origin 
from multiple impurity sources. 

For the given gap, and T = 220K no smaller densities 
than lO^^/cm"^ can be accessed. The optimum njj found 
in Fig.[2lja) translates into n = 1.3-10^^/cm'^, mainly as a 
trade-off between decreasing S"^ (less entropy per carrier) 
and increasing a (larger conductivity with more carriers) . 
In particular, we note that S achieves its maximum for 
smaller concentration of carriers then the figure of merit. 

The thermal conductivity in this range varies very 
slowly with concentration, which is not surprising since 
we fixed the lattice contribution to thermal conductivity 
to a fixed value of kl/Z^ ^ 250W(Km)-i. 

Having chosen parameters to represent FeAs2 , we note 
that the experimentally measured carrier concentration 
in this compound, as inferred from Hall measurements, 
is 5 • 10^''/cm^ in the range of 60-170K,— which is higher 
than the density that optimizes ZT in our model. Thus, 
it seems conceivable that by a deliberate change in the 
impurity concentration or position, an increase in the 
figure of merit of the specimen can be achieved. 

b. lifetime effects. Next we pick the impurity con- 
centration which maximized ZT, and we investigate the 
role of the scattering rate for the figure of merit and the 
Seebeck coefficient in Fig. ^c). The life-time has two 
effects: through the change in chemical potential, and 
directly through the dependence of the response func- 
tions Eq. (fTTj) on scattering rate F. It is this latter effect 
that causes the Seebeck coefficient to vanish at low tem- 
peratures for a sufficiently large scattering amplitude, as 
can be seen in the inset of the top panel of Fig. [2](c) . The 
increase in scattering rate reduces the absolute value of 
both the Seebeck coefficient and the figure of merit, hence 
long lifetimes are preferred in thermoelectric materials. 

The upper limit of the Seebeck coefficient has been 
discussed above. As a function of the scattering rate F, 
figure of merit ZT is limited as well. If lifetimes are long, 
the dependence of the response functions thereof is lin- 
ear, and thus cancels in the dimensionless ratio ZT if 
there are no lattice contributions to the thermal conduc- 
tivity. Therefore, with decreasing F, the figure of merit 
converges towards the purely electronic limit in which 

=0. 

c. position of the impurity level. In Fig. [2jc) we 
show the dependence of the figure of merit on the po- 
sition of the impurity level Eu . We fix the impurity con- 
centration to n^i ~ 2 ■ 10^^/cm^ and the scattering rate 
to F = 5/xeV. It is clear from Fig. [5](c) that a maximum 
ZT is achieved when the donor level is very close to the 
conduction band, which is located at A/2 = lOOmeV. 

This can again be understood from the optimal number 
of carriers : Indeed, in order to reach the ideal electron 
density of 1.3 ■ 10^^/cm^, the chemical potential must 
be rather close to the conduction band. Since at T = 



200K and a gap of A = 0.2eV, the compensating holes 
cannot come from the valence electrons, they have to be 
supplied from the impurity band, hence, for an impurity 
concentration nu = 2 ■ 10^^/cm'^ that is larger than the 
needed number (if completely ionized), also Ed must be 
very close to the conduction band. 

Next we study the sensitivity of the figure of merit 
to the lattice thermal conductivity. With thick line in 
Fig. [2](d) we show the ZT in the absence of lattice ther- 
mal conductivity k l and with a thin line is shown ZT for 
a constant value of kl/Z'^ = 250W/(Km). The figure of 
merit is clearly enhanced when the lattice conductivity is 
reduced, hence the desire for a "phonon glass" (see e.g. 
the review^), i.e. a solid which has a low phonon mean 
free path such as to prevent substantial heat conduction 
by lattice vibration modes. The effect of lattice thermal 
conductivity is most enhanced at low temperature and 
for small impurity concentration. 

In Fig. [2{d) lower panel we plot the optimal carrier 
concentration, which maximizes S or ZT, func- 
tion of temperature. The ZT curve is monotonically 
increasing function of temperature, hence for best per- 
formance at higher temperature we need more impurity 
carriers. To maximize ZT we need larger impurity con- 
centration than we need to maximize 5*. Finally, we also 
display the value of ZT as a function of impurity concen- 
tration for a few representative temperatures. For our 
setup of parameters, the figure of merit is very sharply 
peaked at room temperature around a carrier concentra- 
tion {riD ~ 10-^'^/cm^). 

In conclusion, these model considerations give guid- 
ance as to where to look for promising thermoelectric 
materials. In particular we showed that - unless ver- 
tex corrections or strongly frequency dependent lifetimes 
are of pivotal importance - electronic correlation effects 
are not in the position to enhance the thermopower of a 
gapped system, they can only shift the asymmetry of the 
contributions for electrons and holes. Indeed we found 
the thermopower (of purely electronic origin) to have an 
upper bound that is given by the size of the gap. 



IV. ELECTRONIC STRUCTURE OF FESBz 
AND FEAS2 

FeSb2 and FeAs2 both crystallize in the (regular) 
marcasite structure, have the orthorhombic space group 
Pnnm, and there are two formula units per unit cell (see 
Ref . ISSi ) . The iron ions are surrounded by distorted pnic- 
togen octahcdra, that share corners along the c-axis (see 
e.g. Ref. [33I). In the ligand-field picture, the compounds 
have a Fe3d^ configuration and the 3d-orbitals are split 
into eg and lower lying 2g orbitals. The inequivalence of 
Fe-pnictogen distances causes the 2g to split further into 
two degenerate lower and one higher orbital. In this pic- 
ture, the compounds are in an insulating low spin state 
with the two degenerate 2g orbitals filled 1^ From this 
perspective, metalization of FeSb2 is driven by a tem- 
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(a) FeSba : GGA 



(b) FeSb2 : Hybrid functional 



(c) FcAs2 : GGA 



FIG. 3: Band-structures of FeSb2 (left : GGA, middle : Hybrid functional) and FeAs2 (right : GGA). 



perature induced population of the third 2g orbitali^i^^ 
Previous band-structure calculations, however, suggest a 
more covalent picture^^, in the sense that stabilization 
occurs for d-orbitals that point towards the ligands, i.e. 
in particular lowering the eg orbitals. In FcSb2 this hap- 
pens to the extend that LDA calculation actually yield a 
metallic ground-stat o^^i^^'^^ ' . 

From the perspective of electronic structure methods, 
at zero temperature, the challenge is hence to obtain an 
insulting ground state for FeSb2. We therefore compare 
how three different approaches, GGA, Hybrid function- 
al and the GW approximation perform in this problem. 
In preparation for future work which should include cor- 
relations beyond GW, and to clarify what would be the 
starting Hamiltonian to describe these materials, we ob- 
tain transfer matrix elements and estimate the values of 
the interaction using the constrained random phase ap- 
proximation (cRPA) method. 

d. hand- structures. Our results for the band- 
structure of FcSb2 given by the GGA^ of DFT as 
implemented in the Wien2k package^ are displayed 
in Fig. 13^. We used the atomic positions at room- 
temperature (a=5.83A, b=6.54A, c=3.20A)ii^. Con- 
gruent with previous worksiiiSi2^, the GGA ground- 
state is metallic with small electron pockets half-way be- 
tween the r and Z symmetry points, and corresponding 
hole-pockets at all corners, i?, of the orthorhombic Bril- 
louin zone. Crucial to the understanding of the gap mech- 
anism within more sophisticated techniques (see the GW 
discussion below) is to note the different orbital charac- 
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d^2_y2 d^. 




d.y 


n 


1.27 


1.42 1.42 


1.42 


1.06 



TABLE I: Occupations of the d-orbitals (in the mufBn 
spheres) within the transformed local coordinate system (see 
text for details). The 2g orbitals d^z and dyz are degenerate 
and mainly account for the hole pocket. The electron pocket 



ters of the pockets liiii^ To quantify this, we transform 
the local coordinate system of the d-orbitals into a basis, 
in which the local projection of the d-block of the (GGA) 
Hamiltonian is as diagonal as possible. In this coordi- 
nate system, the x and z axes point (almost) towards the 
antimonide atoms,— and the eg orbitals exhibit the ex- 
pected bonding/anti-bonding splitting. In this basis, the 
electron pocket is mainly of dxy character, and the hole 
pocket is formed by the now degenerate d^z and dyz (2g) 
orbitals. The respective occupations (within the muffin 
spheres) arc shown in Tab. |ll 

The GGA band-structure of FcAs2 (we use a=5.3A, 
b=:5.98A, c=2.88A2^) is shown in Fig. Et. With respect 
to FeSb2, the chemical pressure of the larger As atoms 
is almost isotropic, and the c/a ratio remains virtually 
constant (as a function of external pressure, the ratio 
slightly decreases^). In consequence, and as is apparent 
from the graph, the bands of FeAs2 are much akin to 
those of FeSb2, and could have roughly been obtained 
from a rigid band-shift^. While within GGA the gap 
of FcSb2 was underestimated (no gap at all), the value 
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0.28eV for FeAs2 is just slightly too large with respect 
to the experimental 0.2 — 0.22 eV— i^. On a qualitative 
level, one could thus say that a DFT calculation seems 
to work rather well for FeAs2. Below we will explain why 
we believe this to be a mere coincidence. 

e. Maximally localized Wannier function and cRPA. 
From a conceptual point of view, we find it insight-full 
to note, and compare, the hierarchy of transfer matrix 
elements and the magnitude of local Hubbard interac- 
tions. Starting from a full-potential (FP) LMTO (LDA) 
computation^, we construct maximally localized Wan- 
nier functions for the subsystem consisting of the Fe-3d 
and Sb-5p orbitals, as described e.g. in Ref. lisl - lisl . and 
find that the (largest) nearest neighbor hopping ampli- 
tudes in FeSb2 are tdd ~ 0.2 — 0.3eV, tp^i ^ 0.95eV, and 
tpp ^ 0.7eV. The significance of both the Sb-5p disper- 
sion and the large hybridization was heralded already in 
the band-structure. Fig. and Rcfs. [TtIIisI . as well as in 
the strongly mixed orbitals characters , ^'''^^ see also the 
recent work Ref. |46| . The two transfers tpd, tpp being of 
comparable magnitude, the system thus lies between the 
canonical Anderson model {tpp ^ tpd) and cuprate like 
compounds {tpp <^ tpd). The centers of gravity of the 
d and p bands are separated by Act = 4.2 eV. In FeAs2 
the pd-transfers are larger, as expected from the chemical 
pressure and the covalent /hybridization- like character of 
the gap. We find tdd ~ 0.25 — 0.3eV, tpd ~ l.leV, and 
tpp ~ 0.8eV, and = 5.0 eV. 

To investigate the strength of local interactions, we 
compute the Hubbard U using the constrained RPA^i^ 
technique. Since in FeSb2 and FeAs2 the eigenvalues of 
the pd-subspace are entangled with higher energy bands, 
we employ the scheme presented in Ref. l48l . 

In relation to the rather small energetic shift needed 
to deplete the (GGA) pockets in FcSb2, we find, see 
Tab. ini that the orbital dependence of the Hubbard U 
within the d-orbitals is of notable 5% - an effect hith- 
erto mostly neglected in methods for correlated materials 
that start from a parametrized Hubbard-like Hamilto- 
nian. Before, we alluded to the different orbital char- 
acters of the pockets. Here, we indeed find for the d^z 
and dyz orbital that mainly account for the hole pocket a 
value oi U — 8.6 eV, while for the orbital corresponding 
to the character of the electron pocket the interaction 
is slightly larger with Ud^ = 8.8 eV. This differentia- 
tion in principle favors a charge transfer towards a gap 
opening. Further matrix elements of the interaction are 
Upd ^ 2.5eV, and Upp ^ 4cV. The hierarchy of the inter- 
action strength within the d-shell is already seen in the 
bare, i.e. unscreened. Coulomb interaction and is thus 
linked with the construction of the localized orbitals. 
On general grounds, larger matrix elements of the un- 
screened interaction are indeed expected for orbitals that 
have stronger hybridizations with other orbitals, thus are 
more spatially delocalized^. In analogy with the pres- 
sure dependence of Coulombic interactions in a localized 
basi a^^i^° , the bare interaction of FeAs2 is larger than 
that of FcSb2, in the cubic reference frame we find e.g. 
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7.0 
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dy. 


7.3 


7.2 


7.0 


8.6 


7.3 


d.y 


7.0 


7.5 


7.3 


7.3 


8.8 



TABLE II: constraint RPA values for the Hubbard U (in eV) 
of FeSb2 for the Fe3d-orbitals in the pd-setup of maximally 
localized Wannier functions in the local coordinate system. 



Vdz"^ = 22.7eV and Vdx^-y'i = 22.6eV (the corresponding 
values for FeSb2 are Vdz^ = 22.0eV, Vdx'^-y^ = 21.5eV). 
Moreover, with respect to the antimonide, the Kohn- 
Sham eigenvalues of the arsenide move towards higher en- 
ergies, therewith reducing screening strengths, and caus- 
ing significantly larger values also for the Hubbard U : 



U. 



11.0 eV, U, 



dx^ — y"^ 



10.7 eV. 



/. hybrid functional approach. Previous attempts to 
produce an insulating band-structure for FeSb2 were 
made within the LDA-I-U scheme^, where the param- 
agnetic state (LDA) was found to be stable below a criti- 
cal U=2.6cV with respect to a ferromagnetically ordered 
phase (LDA-f U)ii. Here, we use a hybrid functional ap- 
proach (HYB)i^, and Fig.[3)D displays the resulting band- 
structure, where the B3PW91 functional was used for the 
d-orbitals of the iron atoms^. We also note that, while 
given the freedom, the system docs not develop any mag- 
netic moment within this setup (in LDA-I-U it necessar- 
ily does) . The band-structure features an indirect gap of 
about 0.6eV, i.e. it is by far larger than in experiment. 
This points towards, both, a serious underestimation of 
static correlations within the previously used GGA, and 
the lacking of dynamical effects in the hybrid functional 
approach that will work to reduce the size of the gap. 

g. GW approximation. To investigate the dynami- 
cal effects of electronic correlations in our compounds, 
we applied Hedin's (non-selfconsistent) GW approxima- 
tionSS, which has proven to be quite successful for semi- 
conductors^i^, in its FP-LMTO reahzation^ to both, 
FeSb2 and FeAs2. In Fig. |3^,b wc display, besides the 
FP-LMTO (LDA) Kohn-Sham energies sks, the band- 
structure obtained by taking into account (perturba- 
tively) the energy shifts as provided by the GW self- 
energy 



^KS + Re^{eKs) 



(22) 



with Z-^ = I - i?eE|,, , . In the case of FeSb2 
this indeed opens a charge gap in agreement with ex- 
periment. We note that, consistent with the above hy- 
brid functional calculation, as well as with the discussed 
orbital-dependent interaction strength, the static part of 
the GW self-energy, i.e. setting Z = 1 in Eq. I^^ . yields 
a too large gap of ^ 0.2 cV. Thus, it is the dynamics of 
the self-energy, therewith a true correlation effect, that 
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scales down the gap size with respect to Hartree-hke ap- 
proaches - a situation quite akin to that of correlated 
band insulator a^^'^^ . Indeed, the real part of the diagonal 
matrix-elements of the self-energy are linear in frequency 
over an extended energy range of up to 10 eV. While 
the derivative of these elements are basically orbitally 
independent within the Fe-3d, and Sb-5p orbital subsets, 
respectively, the different hybridizations and also the dif- 
ferent off-diagonal elements yield for the antimonide, in 
the Kohn-Sham basis, a minimal value (eigenvalue of the 
self-energy derivative matrix) of Z « 0.52 for "bands" 
near the Fermi level, and Z « 0.6 — 0.7 for higher ly- 
ing " bands"—. Concomitant with the linear slope of the 
real-part, the imaginary part of the self-energy is basi- 
cally quadratic, but notably asymmetric with respect to 
the Fermi level, see Tab. IIIIl 
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[eV-^] 


FeSba 


FeAs2 
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0.15 


0.08 
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0.02-0.05 


0.02-0.03 



Despite the larger values of the Hubbard J7, the corre- 
lation dynamics is less pronounced in FeAs2, and values 
of Z reach a minimum oi Z k, 0.6 for excitations closest 
to the Fermi level. Also, lifetime effects are both smaller 
in magnitude, and less asymmetric, see Tab. IIIIl While 
for FeSb2 the GW approach has the correct trend with 
respect to experiments, the gap of FeAs2 also slightly in- 
creases from its LDA value, thus departing a bit further 
from the experimental value of ^ 0.2 — 0.22 eV. We note 
that while the gap size within LDA/GGA and GW are 
comparable, the physics is not : In the GW, the size 
of the gap is a result of an almost compensation between 
static (exchange-like) contributions (see hybrids) and the 
dynamical correlations. The induced bandwidth narrow- 
ing within the GW distinguishes its excitations from the 
KS spectrum. 



TABLE III: comparison of the asymmetry of the scattering 
amplitude within the GW approximation. Extraction by fit- 
ting the average d-orbital self-energy (in the Kohn-Sham ba- 
sis) by 3E(|tj| < 5 eV) = -Vu^ . 



V. REALISTIC SEEBECK COEFFICIENTS FOR 
FEAS2, AND FESB2 

In the light of the above considerations for the ther- 
mopowcr of semiconductors, the principle puzzle now is 
why FeSb2, while having a gap that is about 7 times 
smaller than that of FeAs2, has a Seebeck coefficient that 
is (up to) 5 times larger. 

For the calculation of the realistic Seebeck coefficient 
given in Eq. ([T]), we employ the Fermi velocity matrix 
elements of the optics implementation^ of Wien2k, and 
compute the correlation functions according to Eq. ([2]). 



A. FeAs2 

The band-structure underlying the theoretical Seebeck 
coefficient of FeAs2 is the GGA result shown in Fig. [3J:. 
Since the size of the gap is important for the magnitude 
of the thermopower, we scale down its size from its GGA 
value of 0.28 eV to the experimental 0.2 eV. The calcula- 
tion of the transport coefficients uses a small, frequency 
and momentum independent scattering rate/self-energy 
(F ^ 20meV), which therewith practically cancels out in 
the Seebeck coefficient (see the discussion above) and, 
at this point, we do not attempt to introduce effects of 
impurities. 

In Fig. [5] we show our theoretical Seebeck coefficient of 
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FeAs2 as a function of temperature (green dashed curve), 
and compare it to experimental results^^ for the same 
polarization. The agreement is excellent in the intrinsic, 
i.e. not impurity dominated, temperature regime (T > 
12K). 

Also shown is a simple fit, using the formula Eq. (fT7|) 
for the large gap semiconductor. The individual de- 
termination of the parameters 6X and fi (A given by 
experiment) is ambiguous, given the scale of the low 
temperature thermopower (mV/K) with respect to the 
high temperature Heikes limit (which is of the order of 
fcs/e = 86^V/K). In Fig. [5] we show results for the large 
gap model for A/26X ~ n = 85meV (blue dotted curve), 
which is compatible with the constraints \dX\ < 1, and 
H < A/2. 

The decrease of the Sccbeck coefficient at low temper- 
ature can be understood from our model considerations. 
This can both be an effect of the scattering rate as well 
as the presence of impurities, as seen e.g. in Fig. [5Ja,b). 
In the current ab initio case, our limited numerical pre- 
cision (mostly caused by the finite k-mesh) prevents us 
from endeavoring to include these effects. 








T[K] 

FIG. 5: Thermopower of FeAs2. Shown are the theoreti- 
cal Seebeck coefficient for x-polarization, using a constant 
self-energy, and the GGA band-structure, with the gap scis- 
sored to 0.2eV. Experimental results are of Sun et alM' 
(x-polarization). Also displayed is a simple 1/T fit corre- 
sponding to the large gap semi-conductor, Eq. (|17[) . yielding 
A/2SX — fi ~ 86meV. Further indicated is the largest possible 
purely electronic Seebeck coeflicient for FeAs2. 



B. FeSb2 

The situation is entirely different for FeSb2. The 
maximal measured Seebeck coefficient S{T = lOK) is 
— 45mV/K4i If one takes the maximum possible asym- 
metry parameter, = 1, and one assumes the chemical 
potential to be at the most favorable position, i.e. fi = 



—A/2, the charge gap must be larger than A = 0.45eV 
to explain the value of the meassured thermopower in 
terms of our purely electronic model. The experimental 
charge gap, however, is only A w 30meV. We are thus 
led to suspect that the large thermopower of FeSb2 at 
low temperature is not purely of electronic origin. 

A possible scenario, mentioned in the literature, is that 
the very large Seebeck coefficient is mainly caused by a 
substantial phonon drag, i.e. by an electron drift induced 
by a scattering with phonons. While there is no conclu- 
sive evidence that this effect is operational in FeSb2, there 
are several reasons why it is more likely to be present in 
this material than in FeAs2. Since the thermopower is a 
measure for the entropy per carrier, the phonon contri- 
bution to the Seebeck coefficient will be proportional to 
the lattice specific heat times the electron-phonon cou- 
pling constant divided by the electron density. Given 
Debye temperatures of 348 K for FeSb^^^i^s, and 510 K 
for FcAs2^, the specific heat of FeSb2 will be larger than 
that of FcAs2 . The charge carrier concentration at tem- 
peratures where the thermopower is maximal, on the 
other hand, is larger for FeSb2^ : for the best sample 
n 8- 10"/cm^ whereas for FeAs2 n 5 • 10^**/cm^ 

The clcctron-phonon coupling is the least accessible in- 
gredient from the theoretical point of view. Experimen- 
tally there are some insinuations : First of all, the low 
temperature feature seen in the specific heat of FeSb2, 
that has no analogue in the spin response, and is ab- 
sent in the arsenide, could originate from a substantial 
electron-phonon coupling, charge ordering, excitonic or 
polaronic effects from an enhanced coupling to the lat- 
tice. 

Also the nuclear spin-lattice relaxation rate increases 
below 40K,— i.e. in the regime where the thermopower 
starts its huge magnification, whereas an activation law 
decrease (A = 473 K) is found above 50K^ in rough 
accordance with the intrinsic gap. 

Moreover, optical spectroscopy witnesses a large 
change in phonon lifetimes across the metal-insulator 
transition, suggesting an important electron-lattice cou- 
pling^i. Recently, also polarized Raman scattering ex- 
periments gave indications for a notable electron-phonon 
coupling, that is strongly temperature dependent below 
40Ki^ Further, we note that, as expected for substan- 
tial phonon drag contributions to the Seebeck coefficient, 
the magnetothermopower of FeSb2 is very low for those 
samples that exhibit the largest response without mag- 
netic field^. A decrease in the phonon mean free path 
by non-electronic scattering (i.e. in particular by imper- 
fections) is expected to lower the respective effect in the 
thermopower, and indeed the Seebeck coefficient of poly- 
crystalline samples^S and thin filmsH was found to be 
significantly smaller than for single crystals, while hav- 
ing the same high temperature behavior—. Recently, 
also substituted FeSb2_a;As2; was investigatec^^. Inter- 
estingly, it was found that the above mentioned increase 
in the susceptibility, starting at around 50K, is stable 
with respect to the substitution, whereas the shoulder in 
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the resistivity at 10 — 20K is flattened out, and the See- 
beck coefficient decreases. In the phonon-drag picture 
this would, again, be owing to a decrease in the phonon 
mean free path for non-electronic scattering due to the 
presence of the As "impurities" . 
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FIG. 6: Thermopower of FeSb2. Shown are our theoretical to- 
gether with experimental results from Ref.Q for measurements 
along the crystallographic orientations x,y,z, as indicated. We 
can expect reasonable agreement in the temperature range in- 
dicated by the gray area. See text for details. 

Comparing the thermopower of the antimonide and the 
arsenide (see Fig. 1 in Rcf. JJj) one notes that the See- 
beck coefficient of FeAs2 is larger than that of FeSb2 at 
35K and higher. This might indicate - if the phonon 
drag picture holds - that the effective electron phonon 
coupling in FcSb2 has sufficiently decreased (by umklapp 
and phonon-phonon scattering) so that the thermopower 
is now dominated by the electronic degrees of freedom, 
i.e. the larger gap in FeAs2 causes a larger response. Yet, 
we note that optical spectroscopy^! and some transport 
measurements^^i^i see metallic behavior above 70K, or 
already above 40K, respectively, an effect not captured 
by our one-particle approach. Hence we will focus on the 
temperature range from 35K upwards to, at best, 70K. 

Since the GGA Kohn-Sham spectrum is metallic, we 
opted for using the hybrid functional calculation (see 
Fig. [31[b)), albeit with a gap scissored to the experimen- 
tal value A = O.OSeV, to compute the theoretical Seebeck 
coefficient. Moreover, we assume the presence of donor 
impurities at Ed = 9meV, corresponding to an activa- 
tion energy 5 = A/2 — Ejj = 6meV as is seen in the re- 
sistivity in the range of 5-15K)! and we use an impurity 
concentration nu = 10^^/cm'^. This concentration yields 
n'jj + p m 7 ■ lO^^/cm"^ at 20K, in rough accordance with 
the respective hole concentration of 4- 10^ ''/cm'^ found in 
Hall measurements^^. 

We again limit the influence of impurities to their ef- 
fect on the chemical potential.— At high temperature, 
the latter is linear in T as expected, and, using Eq. (PT|) . 



we find an effective mass ratio rjy/rjc = m*^/m* = 0.23 
when using a constant scattering rate, and the very sim- 
ilar rjy/rjc = 0.25 when using the imaginary parts of the 
self-energy from the GW calculation, i.e. the anisotropy 
is mainly propelled by the spectral weight and the Fermi 
velocities, and the GW scattering actually slightly re- 
duces the particlc-hole asymmetry in the current case. 
We further note that the asymmetry is opposite to that 
of FeAs2, where we found rjy/ric ~ 2.5 > 1. 

Thus obtained Seebeck coefficient is displayed in Fig. [SI 
along with experimental results on single crystals^ for 
the three polarizations along the crystallographic axes. 
In the limited range (discussed aboved), starting at 35K, 
and extending towards 70K (indicated by the gray gradi- 
ent in Fig. [6]), we find good agreement with experiments : 
Both, the order of magnitudes, as well as the hierarchy 
of polarizations is captured within our approach. Below 
35K, the single crystal experiment reaches stellar magni- 
tudes of up to -45mV/K— . that we argued to be beyond 
our approach which neglects vertex corrections. Mea- 
surements (not shown) using a polycrystalline sample;^ 
and on films with preponderant < 101 > orientation^ 
display Seebeck coefficients that at low temperatures 
never surpass -500 and -200/iV/K, respectively, while 
having the exact same high temperature behavior, ad- 
vocating a disorder or decoherencc induced lowering of 
the electron drift. At intermediate temperatures, those 
experiments agree qualitatively with both the single crys- 
tal measurements, and our theoretical results. 



VI. CONCLUSIONS 

In conclusion we have considered the problem of ther- 
moelectricity in correlated insulators and semiconduc- 
tors. We developed a simple toy model to study how 
the various many body renormalizations enter the ther- 
moelectric response. We used LDA, hybrid density func- 
tional theory and GW methods to carry out a com- 
parative study of two systems of current experimental 
and theoretical interest FeAs2 and FcSb2. The ratio be- 
tween strength of the Hubbard U and the bandwidth of 
FeAs2 and FeSb2 are comparable and so is the correlation 
strength. In FeAs2 DFT is qualitatively correct, while in 
FeSb2 correlation effects beyond DFT are essential for 
obtaining an insulating ground state, and the one shot 
GW approximation succeeds in that respect. Indeed, us- 
ing this method, we obtained good agreement with the 
experimental values of the gap for both materials. 

The tools developed in this work were sufficient to 
describe the thermoelectric response of FcAs2 quantita- 
tively. This framework is not as successful for the FeSb2 
compound, and in particular it fails to explain the re- 
markably high low temperature thermopower discovered 
by Bentien et alA. Our work implies that the latter 
cannot be understood in the context of local correla- 
tions, and one should focus either on vertex corrections 
to the transport coefficients, or on non local self energy 
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effects characteristic to the proximity to a quantum crit- 
ical point. In this context we notice that within LDA+U 
this material is close to a ferromagnetic instability^. 

An important form of vertex corrections describe the 
phonon drag effect. A framework to estimate quantita- 
tively these effects in conjunction with ab inito methods, 
are currently not available. Above, we mentioned sev- 
eral experimental findings that suggest the presence of 
this mechanism in FeSb2 , providing a strong incentive to 
further development in this vein. 

Future work should include explicit calculations on cor- 
related insulators using LDA-hDMFT to compare with 
the results of the toy model calculations. Furthermore, 
the investigation of vertex corrections on the thermoelec- 
tricity together with the effects of non local self energies 



that go beyond the quasiparticle approximation should 
be considered. 
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of charge neutrality and the particle-hole symmetric point 
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In Ref. IJJ a gap of 0.3eV was extracted from the density 
of states. Yet, the displayed bandstructure is metallic. 
Results (not shown) for low temperature positions^^ do 



not qualitatively differ as far as the following discussion is 
concerned. 

The new unit vectors are x = (—0.627,0.326,0.707), y = 

(0.627, -0.326, 0.707), z = (0.462, 0.887, 0.). 

We note that the former hole pockets on the Brillouin zone 

corners are slightly surpassed as the uppermost valence 

band by momentum regions on the way from the midpoint 

of the xz-face (Y-point), up z- wards towards the edge point 
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